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ABSTRACT 

Using a mass-selected (M* > IO^'Mq) sample of 198 galaxies at < z < 3.0 with HST jmCMOS 
//i6o-band images from the COSMOS survey, we find evidence for the evolution of the pair fraction 
above z ~ 2, an epoch in which massive galaxies are believed to undergo significant structural and mass 
evolution. We observe that the pair fraction of massive galaxies is 0.15 + 0.08 at 1.7 < z < 3.0, where 
galaxy pairs are defined as massive galaxies having a companion of flux ratio from 1:1 to 1:4 within a 
projected separation of 30 kpc. This is slightly lower, but still consistent with the pair fraction measured 
previously in other studies, and the merger fraction predicted in halo-occupation modelling. The redshift 
evolution of the pair fraction is described by a power law Fiz) - (0.07 + 0.04) x (1 -H if -^^^-^. The 
merger rate is consistent with no redshift evolution, however it is difficult to constrain due to the limited 
sample size and the high uncertainties in the merging timescale. Based on the merger rate calculation, we 
estimate that a massive galaxy undergoes on average 1 . 1 + 0.5 major merger from z = 3 to 0. The observed 
merger fraction is sufficient to explain the number density evolution of massive galaxies, but insufficient 
to explain the size evolution. This is a hint that mechanism(s) other than major merging may be required 
to increase the sizes of the massive, compact quiescent galaxies from z ~ 2 to 0. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



Introduction 



The sizes of quiescent massive galaxies at z ~ 2 are 
shown to be on average 3-6 times smaller compared 
to galaxies of similar mass at z - dPaddi et alJl2005l: 
Truiillo et al."2006a"b'; 'Toft et al."2007"; 'Truiillo et al' 



2007; Zirm et al. 20 07; Buitrago et al. 2008; Cimatti et al 
20081; iFranxet alJ r200"8l; lyatTP okkum e tal.l 120081; 



Toft et alJl2009l: Ivan Dok kum et al.i i2010i:iTargett et al. 



201 Ih . High-resolution cosmological simulations con- 



firm the compactness of massi ve galaxies at z ~ 2 com- 
pared to local counterparts dSommer-Larsen & Toft 
2010h . The question is then: what are the physical 



processes that drive the drastic size evolution of mas- 
sive quiescent galaxies between < z < 2? 

At z ~ 2, merging is an important process for the 
evolution of galaxies, in terms of mass /van Dok kum et al" 
,2010) and size: Khochfar & Silk (2006) demonstrate 
with their semi-analytical model that the observed 
redshift-size evolution of elliptical galaxies may be a 
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consequence of the available amou nt of cold gas dur- 



ing the major merger. Furthermore, v an der Wei et al 



1,2009) suggest that major merging is the most impor- 
tant mechanism to produce massive, quiescent galax- 
ies through studying the distribution of the projected 
axial ratio of galaxies from the Sloan Digital Sky 
Survey. On the other hand, there is evidence from 
observation s ( Be zanson et al. 20091) and simulations 
(NaabetaL 2009h that minor mergers are more com- 
mon than major mergers and could be the dominant 
driver for the inferred size evolution. Most of the lu- 
minous red elliptical galaxies at z < 1 are assembled 
through gas- 



s-poor ( i e. dry) rnerging dBell et alJl2004t 
van Dokku J bOOS": 'Bell et al."2006b'). The high frac- 
tion ( ~ 50%,TKi-iek et al. 2006, 2008; IWilliams et al.l 



20091) of massive galaxies at z ~ 2 that are quies- 



cent and have old stellar populations suggests that dry 
mergers may be common since that epoch until z = 0. 
However, it is likely that dry mergers can only account 
for a factor of ~ 2 of growth in size from z ~ 2 to 
dNipoti et al.ll2009h . 

Additionally, gas-rich mergers have been shown 
to drive gas towards the central supermassive black- 
holes and possibly trigger the active galactic nu- 
clei, releasing enough energy to ex pel the gas and 



thereb y quenching star formation dPi Matteo et al 
2005b . The gas inflow can als o enhance star for- 



mation and even fuel starbursts (IBarnes & Hernquist 



19911: iMihos & Hernquist! 119941) . The merging be 



tween two disk-like galaxi es can form an elliptical 
as predicted in simulation s dToomre & Toomrelll972 : 



Barnes & Hernquist 1996"). Although if the merging 



is highly dissipational, a larger degree of rotation and 
therefore flattenirig of the remnant may be expected 
dNaab et all 12006': Robertson et al."2006: 'Wuvt s et al. 



2010l) . Recently, van der Wei et al. (2011) present ev- 



idence for the dominance of such disk-like morpholo- 
gies in quiescent systems at z ~ 2. 

Through studying the abundance of mergers of 
massive galaxies across redshift, we can place con- 
straints on the current evolutionary model of these 
galaxies. Substantial work exists in the literature re- 
garding the merger fraction of galaxies at z < 1.2: 
merger sa mples can ei t her b e constructe d v i a pair 
selection JZepf & K00I 1 19891 Icarlb erg et'aD 1 19941: 



Le Fevre et al 



2000t 



Patton et al.1 



2000; Lin et al 



2004VKartaltepe et al 2007';'Lin et al. 2008;,Bundv et al 



2009; Robaina et al. 2010) or morphological selection 



dLe Fevre et al.l2000l:IConselice et al.l2003l:lLotz et al 
' 2008bl: iMcIntosh et al.l l2008t iHeiderman et al.l 120091 



Jogee et al. l2009h . Kinematic evidence suggests that 
not all irregular morphologies at high re dshift are re- 
lated to mergers (Forster Schreib er et al ]|2009, 2011]). 
Hence, it is not straightforward to identify mergers 
through morphological classification and we focus 
on using pair counts as a probe for merging activ- 
ity in this Paper. It has been challenging to estab- 
lish large samples of pairs of massive galaxies at 
z > 1. High-resolution near-infrared (NIR) imaging 
is required to probe the rest-frame optical emission 
from the stellar populations. Large-area NIR surveys 
have only begun recent l y (e.g., the CANDELS sur- 
vey (IGrogin et al.l 20 1 1 : Koekemoer et al. 20 1 ih . the 



3D-HST survey dvan Dokkum et al.1 12OI iK There 
are only few spectroscopically confirmed merg ers at 



z > 1 .5 dShapiro et al.1 120081; iLaw et al.l 1201 lb . At- 



tempts to constrain the pair fraction at high er red- 



shifts are limited to targetted observations (Bluc k et al 
2009h . The selection criteria of massive galaxies and 



pairs vary across studies, posing a challenge to make a 
uniform comparison of the pair fractions. 

The aforementioned observations compare the ob- 
served pair fraction with the predicted merger frac- 
tion from cosmological simu lations, combined with 
semi-analytical models (e. g .. ISomerville et al.l I2OO8 ; 



Bertone & Consehcd l2009t iHopkins et al.ll2010bh" or 
semi-emp irical models using th e halo occupation dis- 
tribution dHopkins et al.ll2010ah . The potential caveat 
is that these simulations are closely tied to observa- 
tions, often normalized to reproduce the statistical ob- 
servables such as the mass function, the luminosity 
function, and the correlation function of galaxies. 

This Paper uses a sample of 198 massive (M* > 
IO'^Mq) galaxies at < z < 3 with high-resolution 
NIR imaging. The sample is drawn from the Cosmic 
Evolution Survey (COSMOS) where parallel imag- 
ing from the Near-Infrared Camera and Multi-Object 
Spectrograph (NICMOS) onboard the Hubble Space 
Telescope (HST) is available, in order to probe their 
rest-frame optical morphology, off'ering the novel op- 
portunity to derive the pair fraction of a mass-selected 
sample across a wide redshift range with robust pho- 
tometric redshifts (photo-z's) and masses derived from 
spectral energy distribution (SED) fitting. 

The Paper is organized as follows: in ^ we de- 
scribe the photometric catalogue, the quantities de- 
rived and the completeness of the catalogue; the selec- 
tion of galaxy pairs and the correction for projection 
contamination are also discussed. In ^the results of 
the analysis are detailed: we compare our pair frac- 
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tions with other observations and model predictions, 
and estimate the effect of mass completeness on the 
pair fraction. We also explore the impact of merging 
on the growth of the massive galaxy population: the 
merger rates are calculated and the predicted number 
growth is compared with the observed number densi- 
ties of massive galaxies. Implications on our current 
understanding of massive galaxy formation are dis- 
cussed. The conclusions are outlined in ^ 

All magnitudes are quoted in the AB system, un- 
less otherwise stated. A cosmology of //q - 70 km s"' 
Mpc Qm - 0.3 and Qa = 0.7 is adopted throughout 
the Paper. 

2. Data: Catalogue and Selection 



template error function. For sources with /Tveja < 22, 



we mo del the SED in the same way as in iWuvts et al 
(l2007b . in order to estimate the stellar masses. We 



The COSMOS field dScoville et al.ll2007 l) provides 



photometry in 30-1- bands over an area of > 2 deg , in- 
cluding imaging from the Advanced Camera for Sur- 
veys (ACS). The HST/NICMOS Camera 3 (NIC3) 
non-contiguously covers ~ 5 % (332 arcmin^) of the 
field, with 5cr depth of H - 25.6 for point sources. 
The NIC3 imaging is used with the F\6QW filter, and 
the drizzled images have a pixel scale of O.lOl'Vpix 
and a FWHM PSF of ~ 0.25". In our analysis we use 
the NIC3 images from the COSMOS Archiv^Jeduced 
by James Colbert. 

2.1. Input Catalogue 

The analysis of this Paper is based on the pub- 
lic COSMOS 30+ ban d catalogue, combined with the 
i/-band photometry by Gabasch et al.l (120081) . and the 
IRAC p hotometry from s COSMOS. The parent cat- 



alogue dllbert et al. l2009h is selected in the i-band 



(from Subaru Suprime-Cam), where fluxes are mea- 
sured within apertures of 3" in diameter and has a lim- 
iting magnitude of / < 26. The resulting photometric 
catalogue is compiled from all public data in narrow- 
, medium- and broad-bands covering wavelengths in 
UV, optical, NIR and mid-IR, and has a limiting mag- 
nifiide of K < 23.86. 

2.2. Derived quantities 

Photo-z's are derived on all entries using the 
mediu m- and broad-band c atalogue with the code 
EAZY (iBrammer et al.ll2008h . The IRAC fluxes have 
been downweighted by EAZY in the fitting using a 



'http://irsa.ipac.caltech.edu/data/COSMOS/images/nicmos/ 



make use of the BC03 (IBruzual & Charlotll2003l) stel- 
lar populatio n synthesis model with the FAST code 
(iKriek et al. 112009), assuming a Chabrier initial mass 
function (IMF), and fit the SEDs with three diff'erent 
star formation histories: a single stellar population 
without dust, an exponentially declining model with e- 
folding time of 300 Myr and dust attenuation allowed 
to be between Ay - 0-4, and a constant star forma- 
tion model with the same rang e in attenuation. We 
assume solar metallicity and the iCalzetti et al. I (l2000l) 
extinction law. 

2.3. Selection of massive galaxies in pairs 

2.3.1. Parent sample of massive galaxies and mass 
completeness 

We select galaxy pairs by searching for compan- 
ions to the massive galaxies in the NICMOS parallels. 
The parent sample consists of 5,299 massive galaxies 
of M„ > 10" at < z < 3 in the COSMOS field. 
The photo-z's are required to have odds > 0.95 such 
that they have > 95% integrated probability of lying 
within Az = 0.2 of the estimate. The ;t^^-value of the 
SED modelling is required to be less than 10. The 
odds and the^^ criteria reject approximately 55% and 
17% of all the sources at < z < 3.0 in the whole 
COSMOS catalogue, ensuring robustness in the photo- 
z's and masses. 

To estimate the completeness of our adopted mass 
limit of Mi, > 10" Mq from our /-band selected cata- 
logue, we compare the selected gal axies to the .K^-band 
selected FIREWORKS catalogue ('Wuvts et alJl2008h 
for GOODS-Chandra Deep Field South (CDFS), 
which has a deeper limiting magnitude of / = 27 
(3cr). As the completeness is a strong function of red- 
shift, we compare the magnitude distribution of mas- 
sive galaxies in COSMOS and CDFS against redshift 
in Figure [T] Assuming that CDFS is 100% complete 
in selecting massive galaxies, the completeness limit 
of COSMOS are 100%, 75% and 44% for the red- 
shift bins 0-1.7, 1.7-2.3 and 2.3-3.0 respectively. Only 
7% and 14 % of the massive galaxies are rejected by 
the odds and the criteria, so the incompleteness is 
mostly due to the faintness of the massive galaxies in 
the /-band. 

There are 305 massive galaxies in the parent sample 
which have NICMOS i/i6o parallels, but 109 of those 
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are in the edge region of low signal to noise, so there 
are 196 galaxies with usable NICMOS imaging. 

2.3.2. Selection of galaxy pairs 

We run SExtractor Version 2.8.6 dBertin & Arnouts 



19961) on the 10" x 10" NICMOS cutouts (or 35" x35" 
for galaxies at < z < 1 .0), with parameters optimized 
to ensure that sources are deblended properly. The 
isophotal fluxes are used to compute the //-magnitude 
of each source. Due to the large photometric aper- 
ture used in the COSMOS catalogue, 22 cutouts have 
more than one source within the 3"-aperture, where 
the source-confused companions have no separate en- 
try in the catalogue. In these cases the photo-z of the 
companion is assigned to be the same as the primary 
massive galaxy, and the integrated best-fit masses are 
adjusted using the NICMOS //-band flux ratio from 
SExtractor. The final massive galaxy sample con- 
sists of 198 massive (M^, > 10' 'Mq) galaxies at 
< z < 3, where all of them are brighter than the 
depth of the NICMOS imaging described in ^ Note 
that there are two more massive galaxies compared to 
the 196 massive galaxies mentioned in 92.3.11 This is 
because there was source confusion in the photome- 
try of two of the selected massive galaxies, and after 
mass correction there are two massive galaxies on each 
cutout (four massive galaxies in total). 

Galaxy pairs are selected from the massive galaxies 
sample using the following criteria: (1) the massive 
galaxy has one or more companion within a projected 
separation of 30 kpc; and (2) the //if,o-flux ratio of the 
pair is between 1 :4 to 1:1. Imposing these criteria we 
find 40 massive galaxies in pairs (Nobs), in the redshift 
range of < z < 3. Almost all (99%) of the mas- 
sive galaxies are bright enough such that if they have 
1 :4 companions, the companions are brighter than the 
depth of NICMOS (H = 25.6, see g2]i. Only two mas- 
sive galaxies (1%) are fainter than H = 25.6, but they 
are retained in the sample because they have compan- 
ions detected with NICMOS. The number of galaxies 
having companions of flux ratio above 1:2 / 1:3/1:4 is 
20 / 32 / 40 respectively. Examples of the cutouts are 
shown in Figure |2] We note that there are two galaxy 
pairs in which both of the merging galaxies are mas- 
sive (M* > 1O"M0). 

2.3.3. Correcting for chance projection 

Before comparing the pair fraction with model pre- 
dictions and investigating its redshift evolution, it is 



necessary to subtract the contamination from projected 
galaxy pairs at different redshifts. We estimate the 
eff'ect of chance projection by performing a Monte- 
Carlo simulation, assuming that there is no cluster- 
ing in the sources or the massive galaxies. All COS- 
MOS sources are redistributed randomly over the ef- 
fective (unmasked) area of the COSMOS field. The 
198 massive galaxies in the sample are also assigned 
random positions. Using their photo-z's we count, 
within an annulus of 5-kpc to 30-kpc, the number of 
close companions that have magnitudes down to 1:4 
fainter We repeat the redistribution and counting for 
500 realizations. The average of the counts are taken 
as the expected number of galaxies in projected pairs 
(< Npiojected >) for cach redshift bin, and are listed in 
Table [T] It can be seen that approximately half of the 
observed pairs are random projections. We also note 
that the correction is more significant at higher red- 
shift. This is because high-z galaxies are fainter and 
the surface number counts are higher for faint galaxies, 
resulting in a hig her probability of chance projection. 
Law et al. ( 201 lb find a similar correction for chance 
projection (~ 50%) using spectroscopic redshifts avail- 
able for 2874 star-forming galaxies at 1.5 < z < 3.5. 

In each redshift bin, we observe M massive galaxies 
and Nobs of them are in pairs. The fraction of galaxies 
in pairs (fp), or pair fraction for short, is calculated as: 



N 



obs~ 



<N 



projected 



M 



The errors in fp are estimated by the Poisson uncer- 
tainties of Nobs- 

Alternatively, as photo-z's are available for 68% of 
the companions and all of the massive galaxies, we can 
use the photo-z's to reject projected pairs and identify 
the pairs that are physically associated. In practice, 
the pairs are identified by the separation and flux ra- 
tio criteria, and additionally a photo-z critierion: if the 
companion has a separate COSMOS entry with reli- 
able photo-z (odds > 0.95), the 3cr confidence inter- 
vals of the photo-z's must overlap. The number of 
massive galaxies in pairs is given by Np, and the pair 
fraction is given by f, = Np / M. 
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3. Results and Discussions 

3.1. Comparison with previous observations 

In order to examine the redshift evolution of the pair 
fraction, we correct for chance projections in the ob- 
served pairs in COSMOS to get fp in different redshift 
bins through our default approach, as listed in Table [T] 
and plotted on Figure [3] The pair fraction derived us- 
ing the photo-z criterion is remarkably consistent with 
our default approach, except at the highest redshift bin 
in which photo-z's have higher uncertainties and there- 
fore less constraining, but still marginally consistent. 
The Kolmogrov-Smirnov test confirms that the frac- 
tion of massive galaxies in pairs is inconsistent with 
no redshift evolution. We fit the observed fp with a 
power law of the form Fiz) - f (0)(1 + z)'", and find 
the best fit parameters to be F{Q) - 0.07 + 0.04 and 
m = 0.6 + 0.5. 

Robaina et al. I (l2010h (hereafter RIO) use the am- 



plitude of the projected two-point correlation function 
of massive galaxies to estimate the pair fraction of 
galaxies in the COSMOS and COMBO- 17 surveys at 
Z = - 1 .2, requiring galaxy pairs to be separated by 
less than 30 kpc in three-dimensional space, and each 
galaxy to be more massive than 5 x 10'" Mq. Apply- 
ing their mass limit to both galaxies in our pairs, we 
find a fp consistent with their results at z < 1 .2, though 
we note that it is rare to find two galaxies that are both 
massive in a close pair. Our sample of ~ 200 mas- 
sive galaxies yield relatively large uncertainties in the 
pair fraction due to small number counts, compared to 
RlO's sample of ~ 18,000 massive galaxies. The agree- 
ment ensures that our results are com patible with pre- 
viously pubhshed fp below z - \ .2 (IXuetal.ll2004 



Bell et al. 2006a: Mcintosh et al 



200? 



Bundv et al 



20091) that are consistent with RlO's fp. 



Ou r corrected pair fractions agree to that of lLaw et al 
(1201 ih within the uncertainties, though we note that 
their sample are based on star-forming galaxies above 
lO'^Mo. 



Bluck et al.1 (l2009h (hereafter B09) present a study 



of 82 massive galaxies with NICMOS imaging at 1 .7 < 
z < 3 from the GOODS NICMOS survey (GNS), and 
define pairs as any galaxy within 30 kpc and within 
a diff'erence of +1.5 in //leo-magnitude compared to 
the host massive galaxy. They find the pair fraction to 
be 0.19 + 0.07 at 1.7 < z < 2.3 and 0.40 ± 0.10 at 
2.3 < z < 3. The comparison is shown in Figure |3] 
We also compare to the pair fraction from the POWIR 



survey at z ~ 1 from B09. Note that the COSMOS 
catalogue is /-band selected, whereas the GNS targets 
are selected using three different criteria (Distant Red 
Galaxies, Infrared Extrernely R ed Objects and BzK 
galaxies); ( Conselice et al.ll20Tl ). Despite the fact that 
the difference in the selection c ould poten tially bias the 
results, our agreement with the iBluck et al .1 (i2009h re- 
sults is a strong confirmation of high pair fraction at 
z~2. 

3.2. How robust is our pair fraction? 

In our analysis, we quote the pair fraction based 
on the relative fraction of galaxies to a certain depth. 
When the pair fraction is translated to a merger frac- 
tion, it is necessary to account for any systematic bias 
of our massive galaxy sample, i.e. whether the galax- 
ies that we miss due to the limited /-band depth have 
the same pair fraction. We perform a test to estimate 
the eff'ect of the incompleteness on our observed pair 
fraction. In 32.3.11 we demonstrate that the adopted 
mass limit is 75% (44%) complete in selecting mas- 
sive galaxies at 1.7 < z < 2.3 (2.3 < z < 3.0). If 
our sample of 37 massive galaxies is 75% complete at 
1.7 < z < 2.3, we can estimate that in total there are 
~49 massive galaxies, and we miss ~12 of them be- 
cause of their faintness in the /-band. If we assume the 
extreme scenarios, in which all the missed galaxies are 
(not) in pairs, the pair fraction in 1.7 < z < 2.3 would 
then be 0.37 ± 0. 1 1 (0.12 + 0.08). For 2.3 < z < 3.0, a 
similar calculation yields fp of 0.05 + 0.07 and 0.61 + 
0.13 at the limits. It is apparent the conservative lower 
hmits are within the errors of our observed fp. 

The massive galaxies fainter than / = 26 are likely 
to be at the high redshift end, and the faintness can 
be explained by dusty star formation or evolved stellar 
populations. Using the deeper, ^T-band selected CDFS 
catalogue, we find that fp = 0.21+°-2^ at 1.5 < z < 3.0 
for 26 < / < 27, in the fainter regime where the COS- 
MOS catalogue is incomplete. This is consistent with 
our expectation that the missed galaxies would have 
a similar fp as observed in COSMOS for the /-band 
brighter galaxies. We note that the CDFS covers a 
smaller area than COSMOS, and hence statistical er- 
rors in the resulting fp are more severe. 
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3.3. The growth of the massive galaxy population 
through merging 

3.3.1. Merger rate 

We calculate the merger rate as 2l(z) = fp(z)n(z)T^^ , 
where the mer ging timescale (r) is assumed to be 



(A, in units of Mpc ^) is given by: 



r- \^ 

0.4 + 0.2 Gyr dLotz et all l2008al) . and the observed 



co-moving number density of massive galaxies is the 
number of massive galaxies (M) divided by the co- 
moving volume in that redshift range subtended by the 
usable area of 474 NICMOS pointings, i.e. n{z) - M(z) 
/ V'co-moving- The Completeness limits derived in 92.3.11 
are used to correct M. The merger rates are listed in 
Table [T] As an estimate, the uncertainties of /p, n{z) 
and T are approximately 69%, 20% and 50% respec- 
tively. This yields an uncertainty of ~ 88% in the 
derived merger rate. Therefore, the major merger rate, 
unlike the pair fraction, is consistent with no redshift 
evolution within the large range of uncertainties. The 
characteristic time between mergers (F) experienced 
by a galaxy at a given redshift is given by F = t//^, 
and we find the best fit to its redshift evolution to be 
F - 12(1 + z)"'-^. By integrating F over our redshift 
range (see Equation (6) of B09), we estimate that a 
galaxy experiences N,,, = 1.1 + 0.5 major mergers on 
average from z - 3.0 to z = 0, consistent with B09's 
N,„ = 1.7 + 0.5 within the large uncertainties. 

3.3.2. Number density evolution 

The mass function (MF) of galaxies is altered by 
mergers. If the merger fraction at different epochs 
is known, one can translate it into the evolution of 
the massive galaxy population assuming a merging 
timescale (r - 0.4 + 0.2 Gyr in this Paper). We es- 
timate the number of newly created massive galaxies 
using the selected galaxy pairs: for each pair, we cal- 
culate the remnant mass as the sum of the SED masses 
of the galaxies in the pair. In the rare case (4 pairs) 
where the SED mass is not available for the compan- 
ion galaxy because there is no corresponding entry in 
the catalogue, we use the flux ratio and the SED mass 
of the primary massive galaxy to estimate the remnant 
mass. Here we have assumed that the //i6o flux ratio 
corresponds to the mass ratio, and we verify the as- 
sumption by finding consistent remnant masses using 
the flux and the mass ratios for the remaining pairs. 
The number of newly created massive galaxy (Ncreated) 
in each redshift bin is calculated by counting the galax- 
ies that cross the mass limit after merging. The merger- 
induced increment in the co-moving number density 



A = 



^created ^ ^elapsed 



Vco -moving 



X T 



where feiapsed is the time elapsed within the redshift bin. 
Our selected galaxy pairs consists of primary galaxies 
of Mi, > 1O"M0, with companions of flux ratio down 
to 1 :4. Therefore, the remnant mass would be at least 
1.25 X 10"Mo {logM = 11.1). In the case of an equal- 
mass merger, the remnant mass will be 2 x lO^M© 
(logM - 11.3). Normalizing the number density of 
massive galaxies to the observation at z = 2, the re- 
sults are compared with the observed number density 
of massive galaxies above these mass limits, as shown 
in Figure |5] Considering the ~ 0.2 dex uncertainty in 
the number density growth due to counting statistics, 
the slope of the number growth is remarkably consis- 
tent with the observed number density. As the highest 
redshift bin (z > 2.3) is only 44% complete, the pro- 
jected number growth is highly uncertain. The agree- 
ment between our estimated merger-induced number 
density growth and the observed number density sup- 
ports the idea that major mergers are sufficient to ex- 
plain the number density evolution of massive galaxies 
fromz ~ 2.3 to 0. 

One potential caveat of this test is that mergers of 
galaxies less massive than 1O"M0 are not included, 
due to mass incompleteness of the catalogue. An- 
other caveat is the assumption that no new stars are 
formed in the merging, which is only valid for dry 
merging. This can result in an underestimation of 
the number density growth of log(M) > 11.1 galax- 
ies, where equal-mass mergers of two galaxies of 
down to log(M) > 10.8 could contribute to the num- 
ber density. The number density evolution of mas- 
sive galaxies depends on several factors: a merger 
between less massive galaxies can create a massive 
galaxy above the mass limit; on the other hand, if 
two massive galaxies merge, the number of massive 
galaxies would be reduced; the merging timescale is 
closely related to the growth rate of the massive galax- 
ies. The buildup of the massive galaxies can be better 
constrained with larger samples of galaxies down to 
lower masses and higher redshifts, which will be fea- 
sible with the upcoming surveys. The study of num- 
ber density evolu tion is complirnentary to the mass 
density evolution ( IConselice et al ] l2007h and the mass 
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evolu tion of a fixed number de nsity sample across red- 
shift ( van Dokkum et al. 2010l) . in tracing the buildup 
of massive galaxies. Our observed number densi- 



ties show an agreement to van Dokkum et alj (|2010j)'s 
finding that the stellar mass of massive galaxies double 
since z = 2. A precise measurement of the contribution 
of mergers to the number density evolution requires 
accurate determination of the merging timescale, and 
is beyond the scope of this Paper. 

If massive galaxies undergo ~ 1.1 + 0.5 major 
merger between < z < 3 and this is sufficient to 
explain the number density evolution, then this hints 
that major merging can be ruled out as the main mech- 
anism for puffing up the sizes of massive, compact and 
quiescent galaxies from z ~ 2 t o 0, as this size evolu- 
tion requires 2-3 major mergers ( Bezanson et al J2009t 
ToftetalJl2009l) . 



3.3.3. Comparison with models 

To understand how our observations fit into the cur- 
rent understanding of galaxy formation in a cosmolog- 
ical context, we compare our observed pair fraction to 
the expected merger fraction in pair-selected samples 
computed from the "merger rate calculator" (MRC) 
developed by Hopkins et al. (2010a) (hereafter HIO). 
HIO use a halo occupation model to track merger his- 
tory, according to the merger trees constructed from 
the Millennium Simulation (.Fakhour i & Ma 2008). 
The galaxy-galaxy merger rate is determined by con- 
volving the distribution of galaxies in halos with the 
dynamical timesc ale. Assuming a merging timescale 
of 0.35+0. 15 Gyr jLotz et alj2008al) . the merger rate is 
then converted to a merger fraction. Using a simplified 
fitting function, the MRC predicts the merger fraction 
as a function of galaxy mass, gas fraction, redshift 
and mass ratio. We compute the galaxy-galaxy merger 
fraction at < z < 3.0 for galaxies of stellar mass 
between 10" and 10'^ Mq, the range of masses of 
our massive galaxies sample, and of mass ratio down 
to 1:4. The average gas fraction of the pair, defined 
as /gas - Afgas / ( A/gas + ), is a free parameter 
in the model. Direct measurements of the gas mass 
fraction of massive star-forming galaxie s at z ~ 1 and 
Z ~ 2 give 34% and 44% respectively (iTacconi et al] 
2010l) : and 50%-65% for similar systems at z - 1.5 
in another study (Daddi et al. 2010). These systems 
are considered evidence for very gas-rich systems at 
those epochs, hence we select the critical gas fraction 
(/gas) '^o be 20% to differentiate gas-poor ( < /gas < 
/g*Ls ) and gas-rich (/g*, < /gas < 1) mergers. Then we 



multiply the merger fraction, as a function of redshift, 
by two to get the predicted fraction of galaxies in pairs 
to compare with our observations in Figure |4] To in- 
vestigate the importance of dry merging (i.e. nearly 
dissipationless mergers), we overplot the the gas-poor 
and gas-rich pair fraction for comparison. 

Considering that the systematic uncertainties in the 
predicted pair fraction are larger than a factor of two, 
our observed pair fraction is consistent with the predic- 
tion of HlO's model. The number of gas-rich mergers 
are sufficient to explain the number of observed pairs at 
1 < z < 3. Gas-poor mergers are predicted to be more 
frequent than gas-rich mergers below z = 1 (see HIO), 
and are required to explain the observed pair fraction. 
HlO's model predicts 2.1 major mergers per galaxy 
from z = 3 to z = 0, which is almost twice of our re- 
sult and is apparent from the predicted pair fraction in 
FigurelH Our result is also low compared to ~ 1 major 
merger per galaxy at < z < 1.5 for > IO^^^Mq 
(mass limit converted from S alpeter into Ch abrier IMF 
for comparison) predicted bv lDrorv & Alvar ez (200i), 
who estimate the contribution of merging by subtract- 
ing the effect of mass-dependent star formation from 
the galaxy stellar MF. The discrepancy is mostly due 
to the difference in mass limit, confirm ed by a similar 
value o f 1 . 1 major merger predicted by 'Hopkins et alj 
(l2010al) 's MRC if we use Dror y & A lvarez (2008)'s 
mass limit. To reproduce their results of ~ 2 major 
mergers we need to use a merging timescale of 0.25 
G yr, which i s lower, yet still within the uncertainties 
of Lotz et all ( 2008a )'s range of merging timescale for 
pairs having projected separation up to 30 kpc. The ex- 
pected nu mber of major majors fr om D rory & Alvare3 
EqqI and lHopkinsetaP (1201 Oah are based on obser- 
vations of the MF, whereas our result is a direct mea- 
surement of the pair fraction converted into the number 
of major mergers using the merging timescale and the 
number density of massive galaxies. This illustrates 
the need to better constrain the merging timescale, and 
to improve the understanding of how merging alters 
the MF, in order to push merger rate measurements to 
higher accuracy. 

4. Conclusions 

We have quantified the pair fraction of 198 massive 
(M* > 1O"M0) galaxies at < z < 3 from COS- 
MOS with NICMOS parallels. Our findings provide 
a confirmation of the evolution of pair fraction from 
z = 3 to z = 0, in agreement with previous obser- 
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vations ( Bluck et al. 2009b and predictions fro m halo- 
occupation modelling ( Hopkins et al. 2010ah . Gas- 
rich mergers are sufficient to explain the observed pair 
fraction from z = 3 to 1 ; below z - 1 gas-poor mergers 
are also needed. The fraction of massive galaxies ob- 
served to be in pairs is 0.15 + 0.08 from 1.7 < z < 3.0. 
The redshift evolution of the pair fraction is described 
by a power law F(z) = (0.07 + 0.04) x (I + 2)0-6±0-5. 
The merger rate is consistent with no redshift evolu- 
tion, though the uncertainties in pair counts and merg- 
ing timescale restrict the ability to conclusively con- 
strain the merger rate. On average, a massive galaxy 
undergoes ~ 1.1 + 0.5 major merger from z = 3 to 
0, assuming a merging timescale of 0.4 Gyr Using 
the inferred merger fraction, we are able to reproduce 
the observed number density of massive galaxies since 
z ~ 2.3. This implies that major merging can account 
for the number density evolution of the massive galax- 
ies, but other mechanisms such as minor merging may 
be required to explain the size evolution of the mas- 
sive, compact quiescent galaxies at z ~ 2. 
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Table 1 

Pair fraction and merger rate across redshifts 



Redshift range 


No. of massive galaxies (in pairs) 


Expected no. of galaxies 


Pair fraction 


Merger rate %{z) 




M (Nobs) 


in projected pairs <Nprojected > 




(xlO^Gpc-^ Gyr') 


0<z< 1.0 


69 (8) 


2.4 


0.08 ± 0.05 


12.0 


1.0<z< 1.7 


70 (12) 


5.2 


0.10 ±0.06 


7.7 


1.7 <z < 2.3 


37 (12) 


5.9 


0.17 + 0.11 


9.2 


2.3 < z < 3.0 


22 (8) 


5.4 


0.12 + 0.15 


5.6 


1.7 < z < 3.0 


59 (20) 


11.3 


0.15 ± 0.08 
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Figure 1: Mass completeness of the catalogue. Galax- 
ies above our mass limit of M^. > 10" Mq in COS- 
MOS (black) and CDFS (red) are plotted. The gray 
dashed line shows the approximate depth of the COS- 
MOS data. The typical uncertainties in magnitudes and 
photo-z's are overplotted at the bottom left comer. 




Figure 2: The NICMOS /fieo postage stamps of nine 
examples of the selected galaxy pairs. The top row 
shows pairs that were source-confused in the original 
COSMOS catalogue, but are now resolved in our anal- 
ysis with the NICMOS imaging; the bottom two rows 
contain pairs that have individual entries in the cata- 
logue. The IDs and photo-z's of the massive galaxies 
are labelled on the top left and right hand corners of 
each panel. For illustrative purpose, the colour coding 
is scaled logarithmically and the images are smoothed 
by convolving with a Gaussian PSF of FWHM = 2 pix- 
els (0.202"). The angular scale is shown with the 1" 
vertical bar. The white circle overlaid on each map in- 
dicates the 30-kpc search radius around each massive 
galaxy at the centre. 
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■ This work 

This work (projected poirs corrected by photo-z) 
B09 (GNS) 
B09 (POWIR) 
Best-fitting power-low to this work 



0.01 
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Redshift (z) 



2.5 



Figure 3: The redshift evolution of the pair fractions, 
compared to other observations. The black circles de- 
note the fp of our analysis, after statistically correct- 
ing for projection contamination. The black dotted line 
show the best-fitting power law to our /,,, which is of 
the form F{z) = (0.07 + 0.04) x (1 + z)'"'±o.5 -pj^g 
ange diamonds denote our fp, which we use an alterna- 
tive approach to correct for projection contamination 
with the available photo- z's. The gra y squar es and cir- 
cles represent the of iBluck et al.l ( l2009h using the 
GNS and POWIR data. The horizontal bars indicate 
the width of each bin. 
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Figure 5: The redshift evolution of number density of 
massive galaxies. The filled symbols are the observed 
co-moving number density of massive galaxies from 
our sample, with mass limits shown in the legend. The 
lines represent the predicted number growth using the 
observed number density of close pairs, after correcting 
for projected pairs using photo-z. The lines are normal- 
ized to the observed number density at z = 2. 



Figure 4: A plot similar to Figure[3]that compares our 
pair fractions to model predictions. The colour lines 
are the predicted fp for pair-selected samples f rom the 
merger rate calculator of lHopkins et al. I ( 12010 j) assum- 
ing different sets of gas fraction. 
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